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We investigate the consequences for geometrically frustrated antiferromagnets of weak disorder in 
the strength of exchange interactions. Taking as a model the classical Heisenberg antiferromagnet 
with nearest neighbour exchange on the pyrochlore lattice, we examine low-temperature behaviour. 
We show that random exchange generates long-range effective interactions within the extensively 
degenerate ground states of the clean system. Using Monte Carlo simulations, we find a spin glass 
transition at a temperature set by the disorder strength. Disorder of this type, which is generated 
by random strains in the presence of magnetoelastic coupling, may account for the spin freezing 
observed in many geometrically frustrated magnets. 
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Geometrically frustrated magnets - materials in which 
magnetic ions form a lattice consisting of frustrated units 
such as triangles and tetrahedra - characteristically re- 
main in the paramagnetic phase even at temperatures 
that are small on the scale set by exchange interactions. 
Nevertheless, at sufficiently low temperature most exam- 
ples exhibit spin freezing, with a large value (in the range 
from 10 to 150) for the ratio |Ocw|/^F of the magnitude 
of the Curie- Weiss constant to the freezing temperature 
as an identifying feature [l[. Indications of freezing in- 
clude: a difference between field-cooled and zero-field- 
cooled susceptibilities 0, 0, 0, H, @, 0] ; a suppression of 
inelastic magnetic neutron scattering @, H, H, @] ; an d m 
some cases a divergence in the non-linear susceptibility 
on approaching Tp, as in conventional spin glasses [1 01 ] - 
Some well-studied systems are SrCr8-a;Ga4+a;Oi9 (with a 
layered structure consisting of slabs from the pyrochlore 



lattice) @, i, 
mum jarosite 



the kagome antiferromagnet hydro- 



t nyd.ro- 

|, and the pyrochlore Y2M02O7 @, 0]. 

The origin of this spin freezing has long been puzzling. 
On the theoretical side, it is established that freezing is 
absent from some simple models without disorder, includ- 
ing the classical Heisenberg antiferromagnet with near- 
est neighbour interactions on the pyrochlore lattice 11 1 . 
While a more realistic treatment should take account of 
disorder and of various residual interactions, in some of 
these materials there appears to be little structural dis- 
order Q, and in others Q Tp is rather insensitive to the 
identified form of disorder, dilution at magnetic sites. 

In this context, recent experiments which show the 
importance of random strains in two pyrochlore mate- 
rials are particularly interesting, since via magnetoelas- 
tic coupling such strains will lead to local variations 
in the strength of exchange interactions. One mate- 
rial is Y2M02O7, in which disorder in Mo— Mo dis- 
tances has been revealed using XAFS (T^. The other 
is Zni_ a; Cd a ;Cr204. Unusually for a geometrically frus- 
trated magnet, in pure form (x=0) this has a first-order 
transition to a low temperature phase in which magnetic 
degeneracy is lifted by a lattice distortion and there is 



Neel order [13J. This phase has a striking sensitivity to 
substitution of Cd for Zn: x—0.03 is sufficient to suppress 
Neel order completely, with spin freezing taking its place 
(l3| . The effect is argued 14 1 to arise from strains around 



Cd sites, generated because Cd 2+ has a larger ionic ra- 
dius than Zn 2+ . In remarkable contrast, Neel order is 
much less sensitive to magnetic dilution by substitution 
of Gd for Cr: it survives with up to 25% of magnetic ions 
removed flij ]. 

Taking the background outlined above as motivation, 
our aim in this paper is to discuss the effects of weak 
exchange randomness in geometrically frustrated antifer- 
romagnets. We focus on the classical Heisenberg model 
with nearest neighbour interactions on the pyrochlore lat- 
tice, because of the large body of experimental work on 
pyrochlore antiferromagnets, and because its behaviour 
without disorder is well understood 11[. In particular, 
the disorder-free model is known to have an extensively 
degenerate, connected ground state manifold [llj, and 
a dipolar form for s pin correlations in the limit of low 
temperature T 
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disorder one expects two regimes according to its am- 
plitude. Characterising interactions by their average J 
and a magnitude A of fluctuations, for strong disorder 
(A > J) the model is simply a conventional example of 
a spin glass. Our interest lies instead with the weak dis- 
order limit (A <C J) in which exchange randomness acts 
as a perturbation lifting the ground-state degeneracy of 
the clean system. In the following we show that projec- 
tion of fluctuations in nearest-neighbour exchange inter- 
actions into the ground state manifold of the clean system 
generates long-range effective interactions. In addition, 
using Monte Carlo simulations with parallel tempering, 
we show for A <C J that the model has a transition at 
a temperature T-p cx A. As the transition is approached 
from above, the spin glass susceptibility diverges, and for 
T < Tp there is long-range spin glass order. 
We start from the Hamiltonian 



n 



ij 



(1) 
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elegant way of building in that constraint is to parame- 
terise the ground states in terms of a gauge field, as set 



of weak exchange disorder can be introduced into that 
formulation. The essence of the gauge field parameteri- 
sation for the model without disorder can be summarised 
as follows. A set of vector fields B a (r) is introduced to 
represent spin configurations, with one field for each spin 
component. The ma pping b etween a spin configuration 
and vector fields (see 15|,LL6j for details) is made in such a 
way that the condition for a configuration to be a ground 
state translates into the condition V ■ B(r) = 0. After 
coarse-graining, the fluxes B a (r) are treated as continu- 
ous, divergence-free fields, with a statistical weight e~ s ° 
(before normalisation) given by 



Do 



d 3 rV 



B Q (r) 



(3) 



in which classical spins Si are three-component unit vec- 
tors at the sites i of a pyrochlore lattice, and exchange 

interactions Jy are non-zero only between neighbouring out in Refs. 15|, . We next examine how the effects 
pairs of sites. We begin with a qualitative discussion of 
the special features of this model at weak disorder. 

As a first step, consider a single tetrahedron taken 
from this lattice, with spins Si . . . S4 at the vertices. Its 
ground states are the configurations for which J2i Si = 0. 
With all Jij equal, spin stiffness is zero in this toy prob- 
lem in the sense that, within the set of ground states, 
the orientations of a pair of spins can be chosen arbi- 
trarily. The consequences of small variations in Jy with 
amplitude A have been set out in Ref. [l7j]: generically, 
a unique ground state is selected (up to global spin ro- 
tations) in which all four spins are collinear and the spin 
pairs linked by the strongest interactions are arranged an- 
tiparallel. Variations in Jy hence induce a ground state 
stiffness, since changes in the relative orientation of a pair 
of spins cost an energy 0(A). 

Moving to the full, pyrochlore lattice problem, we next 
argue that weak exchange randomness generates long- 
range effective couplings. The logic is as follows. Without 
disorder Jy can be block-diagonalised by Fourier trans- 
form. Its spectrum has four branches in the Brillouin 
zone. The lowest two branches are degenerate and in- 
dependent of wavevector, mirroring the ground state de- 
generacy of the model. In the limit A <C J it is natural 
to project the matrix Jy onto this degenerate subspace. 
The matrix elements Py of the projection operator have 
a dipolar form 18): they decrease as |rj — r?| -3 for large 
|rj — Tj\. Hence so does the projected interaction ma- 
trix. In the context of conventional spin glasses, dipo- 
lar interactions have been shown not to be sufficiently 
long ranged to generate different critical behaviour from 
that with short-range interactions [lj|; it is not clear at 
present whether this conclusion carries over to the prob- 
lem we are concerned with. 

The interaction matrix and, in the limit A « J. its 
projected version, enter directly into a calculation in 
which the fixed spin length | | 2 = 1 is treated within 
the spherical approximation X)i|Si| 2 = N s (where the 
sum is over N s spins in the lattice). In the absence of 
disorder such an approach is exact for a model in which 
the number of spin components n — > 00 [20], and give s 
an excellent treatment of low T correlations at all n [l5| . 
Applying it and the replica method, with a Gaussian dis- 
tribution for Jij of variance A 2 , we find [2l[ for A <C J 
that there is a spin glass transition at a temperature in- 
dependent of J and proportional to A: 



where the stiffness re is determined by microscopic de- 
tails of the model. From the example of a single tetra- 
hedron (and the details of the mapping between spins 
and vector fields), we know that the effect of small vari- 
ations in is to favour a specific axis for flux. This 
axis turns out to be one of the cubic crystal axes: the 
particular one selected depends on the values of Jy and 
varies randomly from one tetrahedron to another. We 
therefore propose an effective theory for a geometrically 
frustrated antiferromagnet with weak exhange disorder, 
taking S e fi = So + Sdis + Smt with 



0A J d 3 r^[B a (r)-n(r)f 



and 



Sdis — 



S int =(3u I d 3 r[V|B a (r)| 2 ]' 
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(5) 
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A difficulty in proceeding further with a conventional 
replica treatment lies in restricting spin configurations 
to the manifold of ground states for the model without 
disorder, as is necessary for A, T <C J- A natural and 



Here, the preferred local axis for flux is defined by the 
random field n(r), which has average [rij(r)] av = and 
variance [ni(r)rij(r')] av = <5y<5(r — r'). Inverse temper- 
ature is denoted by f3: Sdis dominates over So at low 
temperature, since the first is an energetic contribution 
while the second is entropic. Microscopies imply an up- 
per bound to |B a (r)|, imposed here by Si nt with phe- 
nomenological coefficient u. Note that S e ff, like 7i is 
invariant under global spin rotations, which translate to 
rotations between the different B a (r) . In the langauge of 
this effective theory, spin freezing is condensation of flux 
into a specific arrangement favoured by disorder. 

To search for such freezing, we turn to Monte Carlo 
simulations for the model of Eq. Our work builds 
on the initial investigation of Ref. [22f . but is much 
more detailed. In common with a recent study of the 
spin glass transition in the three-dimensional Edwards- 
Anderson Heisenberg model [23[ , we use parallel temper- 
ing 24 1 to reach equilibrium at low temperature. Some 
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FIG. 1: Evolution of C(r) with Monte Carlo time (note log- 
arithmic scale) for T/J = 0.01, A/J = 0.1, L = 7, and max- 
imum r; upper and lower curves from different initial states, 
as described in text. 



details are as follows. We take Jy independently and 
uniformly distributed on the interval [J — A, J + A] with 
0.05 < A/ J < 0.2, and study the temperature range 
10~ 2 < T/J < 1. System sizes, specified by the lin- 
ear dimension L in units of the lattice constant and by 
the number of spins N s = 4L 3 , are 2 < L < 7 and 
32 < N s < 1372. Run lengths vary from 5 x 10 3 Monte 
Carlo steps per spin (MCS) for L = 2 to 2 x 10 5 for 
L = 7. Results are averaged over a number of disorder 
realisations varying from 10 3 for L = 2 to 200 for L = 7. 

We present data for three quantities: the heat capacity 
per spin C v , the spin glass correlation function C(r), and 
the spin glass susceptibility %. The last two arc defined 
in terms of the behaviour of two copies of a system with 
identical disorder. Labelling the copies with I = 1,2, 
denoting thermal averages in each copy by (...}/ and a 
disorder average by [. . . ] al) , we have 



and 



C(r) = [(S(0).S(r)) 1 (S(0).S(r)) 2 ] 



(6) 



(7) 



First, as a test for equilibration, we examine C(r) as a 
function of simulation time, comparing initial conditions 
for which C(r) = 1 (both copies initially in the same 
Neel ground state of the diorder-free model) with ones 
for which C(r) = (each copy initially an independent 
random spin configuration). Results are shown in Fig. [I] 
for the most demanding case (maximum L and r, mini- 
mum T). On this basis, for L = 7 we collect data after 
discarding an initial 1 x 10 5 MCS. Equilibration is much 
more rapid at smaller L or larger T. 

The heat capacity, shown in Fig. [2j varies smoothly 
with T, as expected in the absence of a Neel ordering 
transition. The consequences of exchange randomness 




FIG. 2: C v (in units of k B ) vs T/J for A/J = 0.1 and (from 
top to bottom) L—6, 5, 4 and 3. 

are revealed in the limiting value of C v at small T: with- 
out disorder this is [ll| 3&b/4 for large L, reflecting the 
zero modes (1/4 of all degrees of freedom) of the ground 
states, but with disorder it rises to /cb, because all zero 
modes are removed, except for the three arising from 
global spin rotations. 

Next, we present the central result of our simulations, 
the behaviour of the spin glass correlation function C(r), 
shown in Fig. [3l We believe this provides clear evi- 
dence that C(r) is non-zero at large r below a transition 
temperature T F . From simple inspection of this figure 
0.04 < T F < 0.02 at A/J = 0.1. In an effort to determine 
T F more precisely, and to illustrate finite-size effects, we 
turn to the susceptibility x, shown in Fig. [4] The rapid 
increase in \ with L at low T is clear. Close to Tp one 
expects the finite-size scaling behaviour 



X (T,L) = Li/»f(L 1 /»t) 



(8) 



where t — (T — T c )/T c , and v and 7 are the standard 
critical exponents for the correlation length and suscep- 
tibility. Scaling collapse of the data is shown in the inset 
to Fig.H with T F = 0.023, v = 1 and 7 = 1.45. Uncer- 
tainties in these parameters are hard to quantify because 
finite size effects are large for L ~ 2 and L — 3 (N s = 32 
and N s — 108); omitting these sizes, for A/J = 0.1 col- 
lapse is obtained with 0.020 < T F < 0.32, 0.9 < v < 1.2 
and 1 < 7 < 1.6. 

Finally, we examine the dependence of Tp on A. On 
dimensional grounds, one has T F /A = g(A/ J), and from 
our discussion of spin stiffness in a single tetrahedron we 
expect g(x) to be finite in the weak disorder limit x — > 0. 
We have evaluated C(r) as a function of T for A/ J = 0.2, 
0.1, 0.75 and 0.05, considering only L = 5 (N s = 500) 
because of limitations on computational resources. As 
shown in Fig.[5j the dependence of data for C(r) at fixed 
r on T and A can be reduced to a single scaled variable 
T/T (A), and T (A) oc A for small A. 

In summary, we have shown that weak exchange ran- 
domness in the classical Heisenberg antiferromagnet on 
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FIG. 3: C(r) vs r for L = 7 and A/J = 0.1 at (from top to 
bottom) T/J = 0.01, 0.02, 0.04, and 0.1. Inset: dependence 
on system size at T/J — 0.01, L = 5, 6, 7. 
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FIG. 4: x vs T/J for A = 0.1J and system sizes from L — 2 
to L = 7. Inset: scaling collapse. 

the pyrochlore lattice generates long-range effective in- 
teractions, and that these are responsible for a spin glass 
transition at a temperature set by the disorder strength. 
We suggest that this may account for spin freezing ob- 
served in many geometrically frustrated magnets. 
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FIG. 5: Dependence on disorder strength A: C(r) for r = 4 
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guide to eye. 
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